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Abstract. Various empirical and theoretical studies indicate that cumulative network traffic is a Gaussian 
process. However, depending on whether the intensity at which sessions arc initiated is large or small relative 
to the session duration tail, Mikosch et al (2002) and Kaj and Taqqu (2008) have shown that traffic at large 
time scales can be approximated by cither fractional Brownian motion (fBm) or stable Levy motion. We 
study distributional properties of cumulative traffic that consists of a finite number of independent streams 
and give an explanation of why Gaussian examples abound in practice but not stable Levy motion. We offer 
an explanation of how much vertical aggregation is needed for the Gaussian approximation to hold. Our 
results are expressed as limit theorems for a sequence of cumulative traffic processes whose session initiation 
intensities satisfy growth rates similar to those used in Mikosch ct al (2002). 



1. Introduction 

Collection of data network measurements often uses an algorithm for clustering packets with the same 
source and destination IP addresses. Various criteria for grouping packets yield different entities, e.g. con- 
nections, flows (or unidirectional connections), end-to-end streams, etc. (See e.g. Sarvotham et al, 2005, 
Section 4). These high-order constructs of packet clusters are sometimes termed sessions. For now, think of 
a session as a user downloading a file, streaming media, or accessing websites. For each session, summary 
measurements are computed for the size (the number of bytes transmitted in a session) , the duration of the 
session and the average transfer rate. Data sets of these summaries show some distinctive properties, such 
as heavy tails for session size and duration (Arlittson and Williamson, 1996; Crovella and Bestavros, 1997; 
Willinger et al, 1998) and sometimes rate (Maulik et al, 2002; Resnick, 2003). 

Typically, a time resolution or granularity is selected or imposed. Typical resolutions are 1, 10 or 100 
milliseconds, 1 second, 1 minute, 1 hour, etc. Once a resolution is fixed, the number of bytes or number of 
packets per unit time can be recorded and cumulative network loads over stationary time intervals computed. 
These cumulative loads have been studied from empirical and theoretical perspectives with the objectives of 
satisfying performance criterion and offering adequate bandwidth provisioning (van de Meent and Mandjes, 
2005) or predicting properties of congestion events (Jin et al, 2007). 

Conventional wisdom based on empirical studies claims that a heavily loaded network link subject to 
aggregation over many users should see Gaussian traffic. This wisdom is considered a network invariant. 
Influential examples based on the Bellcore measurements (Leland et al, 1994) suggest that horizontal aggre- 
gation, that is, working with a single on/off stream at sufficiently large time scale justifies Gaussian modeling. 
See also Kurtz (1996) and Willinger et al (1997). 

However, mathematically it is known that with heavy tailed session durations, cumulative load at large 
time scales can be approximated by either fractional Brownian motion (fBm) or stable Levy motion, de- 
pending on whether the intensity at which sessions are initiated is large or small relative to the size of the 
duration tails. See Mikosch et al (2002); Kaj and Taqqu (2008); Taqqu ct al (1997). The stable approxi- 
mation has not been observed empirically (Gucrin et al, 2003) and use of Gaussian cumulative loads has 
become dominant (Kilpi and Norros, 2002; Sarvotham et al, 2002; Jain and Dovrolis, 2005). 

But why should traffic be Gaussian? According to the empirical study van de Meent et al (2006), in addi- 
tion to horizontal aggregation, the superposition of independent traffic streams, that is, vertical aggregation, 
can justify a Gaussian model and, in fact, the number of traffic streams need not be large to make cumulative 
loads approximately Gaussian. 
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In this paper we 

• study the distribution of the cumulative load in the presence of a finite number of independent traffic 
streams; 

• give an explanation for why Gaussian examples abound in practice but not stable ones; 

• answer how much vertical aggregation is needed to justify the use of fBm. 

Our findings suggest that cumulative load for aggregate traffic can be approximated by fBm at large time 
scales provided the initiation intensity of at least one of the traffic components is large. Network traffic in 
the wild has several distinct constituents and we claim that in practice there is one or more components 
with dominant large initiation intensities. For example, this should be the case with web traffic using port 
80 and this suggests why Gaussian traffic should be pervasive (van dc Mecnt ct al, 2006). 

Before discussing mathematical details, we illustrate the phenomena of interest with a motivating example 
of a network trace captured at Cornell University main campus servers during 55 days between November 2, 
2009, and January 15, 2010. Cornell's data set is a collection of netflow records, where all non-IP traffic has 
been discarded and only TCP and UPD traffic is present in the trace. A netflow is a collection of packets 
with the same source and destination IP addresses, source and destination ports, protocol, ingress interface 
and IP type of service (Cisco Systems, Inc., 2007). In our data, TCP traffic accounts for nearly 90% of the 
bytes, and over 80% of the total number of netflows, mostly port 80 (http traffic) netflows. We have taken 
the part of the trace corresponding to both outgoing and incoming traffic between I and 5 p.m. local time, 
adding up to 220 hours of traffic. The anonymization procedure used on the data obliterated the distinction 
between outgoing and incoming flows. 

We analyze the distribution of A^ TCP ^ and A^ UDP \ namely the cumulative load generated by TCP and 
UDP bytes, respectively. For this purpose, we separate the trace into TCP and UDP netflows and for 
k = I , . . . , 220 we count 

(TCP) 

A k := total number of TCP bytes captured in the fcth hour, 
D := total number of UDP bytes captured in the fcth hour. 

Due to the dates and times of collection, these counts exhibit both a trend and a daily seasonality. Here we 
detrend and remove daily seasonality (see e.g. Brockwell and Davis, 1991, Section 1.4), but our conclusions 
are the same without this massage. 

Figure 1 shows Gaussian QQ plots for A^ TCP ^ {left) and A^ UDP ^ {right). A straight line fit is evident for 
the TCP cumulative input. However, the UDP counterpart shows a significant departure from the straight 
line. Using the p- values of the Anderson-Darling two-sided test also shows no evidence against the normality 
of A^ ^ {p = 0.1369), but strong evidence against a Gaussian model for A^ UDP ^ {p = 9.8 x 10~ 16 ). 

We also check whether A^ UDP "> is a heavy-tailed random variable, in the sense of its distribution tail 
being regularly varying with tail index a (de Haan and Ferrcira, 2006; Resnick, 2007). For instance, Figure 
2 left shows a stable regime in the Hill plot of a (for Hill plots, see e.g. Hill, 1975; de Haan and Resnick, 
1998; Resnick, 2007). Additionally, in Figure 2 right we present the exponential QQ plot of \og{A^ UDP ^ ) 
with a straight line fit through the biggest 55 observations. This shows no evidence against approximating 
the distribution of thresholded values of A^ UDP ^ by a Pareto (Recall that the logarithm of Pareto random 
variable is exponential; sec, for example, Mc Neil et al (2005); Resnick (2007); Coles (2001).) 

If we consider the aggregated cumulative load, A^ TCP ^ + A^ UDP \ the normal QQ plot in Figure 3 exhibits 
a straight line fit and the Anderson-Darling test p— value is 0.2117, showing no evidence to reject normality. 
Without accounting for centering and scaling, this result is rather counterintuitive due to the nature of the 
individual tails of A^ TCP ^ and A<- UDP K 

Our explanation to the above phenomenon starts by modeling the quantity of data in windows of length 
T in Section 2. Analogously to the slow and fast growths of Mikosch et al (2002), we define two different 
scenarios for the aggregated traffic. A third scenario is defined similarly to the boundary case considered in 
Kaj and Taqqu (2008). In Section 3 we obtain approximations and provide clarification of the asymptotic 
behavior at large time scales. We let T — > oo and see what limits exist for the aggregated cumulative load. 
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FIGURE 2. Plots for the UDP cumulative input. Left: Hill plot of tail index with 95% 
confidence interval. Right: Exponential QQ plot of the log data. 



In Section 4 we study extensions to our model and finally Section 5 contains some technical results used to 
prove our main theorems. 

1.1. Notation. In order to simplify the later presentation, we introduce and collect some notation. Refer- 
ences are provided for further reading. 
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Figure 3. Normal QQ plot of the aggregated cumulative input. 
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The right tail of the distribution function F, i.e. F = 1 — F. 
The left continuous inverse of the distribution function F, 
i.e. F*~(y) = inf{a; : F(x) > y}. 
lim^oo f\{x)/f2(x) = 1. 

Convergence of finite dimensional distributions. 

Vague convergence of measures. See e.g. Kallenberg (1984); Resnick (1987). 
The space of nonnegative Radon measures on (0, oo]. 
The class of regularly varying functions with index 7. 

See e.g. Bingham et al (1987); de Haan and Ferreira (2006); Resnick (1986). 
A Poisson random measure £ with mean measure F£. 

o 

A compensated Poisson random measure with mean measure E£, i.e. £ = £ — F£. 



N, 



j U h+^du ■ hs(dr)) on M. x (0, oo) 2 , where hs is a measure on (0, 00) 



such that h§[-, 00) € RV-g- If h$(dr) = Sr we may simply write N?f 



7,(5' 



A^ 



- ^ PMM(ds ■ ju'^+^du ■ 8r-( s+1 Ur) on K x (0, oo) 2 



A 7— stable random measure with control measure m(dv) and stable index 1 < 7 < 2. 
For £ = PKM(m(dv)w-( p+1 ^dw), we can write 

M„ m (A) £ ((-cob 2?) f^) _1/T / A C=o^(^^)- 
See e.g. Samorodnitsky and Taqqu (1994, Chapter 3). 

A 7— stable Levy motion totally skewed to the right with stable index 1 < 7 < 2. 

In general, we can write A 7 (i) = f (— cos ^) ^ 
See e.g. Samorodnitsky and Taqqu (1994, Chapter 3). 
The standard fractional Brownian motion with Hurst exponent H 



)i/7 
jo°° 1 {o<v<t}M p>m (dv). 
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2. Model Description and Basic Assumptions 

Consider a network that has an infinite number of nodes. At certain times, a node begins a transmission 
session at a random rate that is fixed throughout the session. Suppose network traffic consists of p distinct 
types which we call streams. In practice, such a division of network traffic arises naturally; e.g. traffic can 
be segmented by application type (web, email, streaming media, file-sharing applications, etc.), by protocol 
(TCP, UDP, IMTP, etc.), and even by users. We suppose the p streams are independent and that each 
follows an M/G/oo input model. The overall load is obtained by aggregating over the p streams. Thus, the 
basic assumptions are as follows: 

• Sessions corresponding to the jth stream are initiated at homogenous Poisson time points 
{r^, — oo < k < 00} with arrival intensity X^ > 0. These points are labeled so that < < t[ 3) 
whence {— Iq ,T± \ (rjf^ — rjj? , k ^ 0)} are iid exponential with parameter Thus, we have: 



= FRM(X U) ds). 



We assume that these PRMs are independent. 

• All the sessions in the network transmit data at positive random rates that are iid with common 
distribution Fr. Let {iil } be the rate of the fcth session of the jth stream. Assume that either 
F R E RV- aR , 1 < a R < 2, or E^R^) 2 ] < 00. In either case, define fj, R := ER^> . 

• Sessions in the jth stream have positive durations {-Dt 7 }, j = 1, ...,p, that are iid F%\ with 
F$ € RV_ o) , 1 < ay' < 2, and p£> := ED[ 3) . In general, not all the aQs are equal. 

• We also assume mutually independent durations across streams, and that durations and rates are 
independent. 

There is empirical evidence justifying the choices of ap's and a^: See e.g. Cunha et al (1995); Willinger et al 
(1995); Lcland et al (1994); Resnick (2003); Lopcz-Oliveros and Rcsnick (2009). For now, we adopt a 
network-centric approach by assuming the rate of communication entirely depends on the state and speed of 
the network. Studies supporting this assumption include Shakkottai ct al (2005) and Kortcbi et al (2005). 
We will need 

A = J>Cfl (2.1) 
3=1 

F D :=52(\W/\)Fg, (2.2) 

and the quantile functions 

= (i/4 j) rw - (p^ra i/t), (2.3) 

b D (t) = (l/F D r(t) - F£(l - l/t), (2.4) 
b R (t) = (l/F R r(t)=F R -(l-l/t). (2.5) 

Notice that Fd is the mixture model of the durations of the p streams, with weights A^-* /A, j = 1, . . . ,p. 
In fact, Fd is the distribution of the duration of the sessions of the aggregated stream, and X^'/X is the 
proportion of the traffic that consists of sessions from the jth stream. We return to this interpretation later. 

Now consider (s, u, r) as a generic Poisson point representing a session that starts at time s, has duration 
u and rate r. By augmentation, the counting function of the session descriptors (rj^, Dj?\ R^) of the jth 
stream on K x [0, oo) 2 is 

NU) ■=J2^KD^\R^)= im ( xU)dsF D(du)FR(dr)), j = l,...,p. (2.6) 

k k k k 
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By independence, the counting function of the session descriptors of the aggregated stream is 

v I p 

N :=J2 Nii) = PMM Ads^(A (j) / 'X)F { J ) {du)F R (dr 
i=i V 3=1 

= ¥WM(XdsF D (du)F R (dr)). (2.7) 

Thus, the mean measures of the N^' and N are given by 

EN u) (ds,du,dr) := A (j) dsF^ (du)F R (dr), j = 1, . . . ,p, 
EN(ds,du,dr) := XdsF D (du)F R (dr). 

In addition, let 

rt ru 

L t (s, u) = |[0, f] n [s, s + u}\ = / l [s , s +u]{y)dy = l [0 ,t](y + s)dy, (2.8) 

Jo Jo 

be the length of the subinterval of [0, t] during which the session (s, u, r) transmits data. In Lemma 5.4, we 
summarize several required properties of L t (s, u). 
For each j, define 

A^'(t) := cumulative input in [0,t] from the jth stream 



oo poo poo 



rL t {s,u)N ( - : > ) (ds,du,dr), (2.9) 

J-ooJO JO 

and similarly 

A(t) := cumulative input in [0,t] from the aggregated stream 

rL t (s,u)N(ds,du,dr). (2.10) 



oo poo poo 



-oo JO 

(Kaj and Taqqu (2008) showed that these integrals are well defined using Campbell's theorem (Kingman, 
1993, Section 3.2).) Also, 

EAW{t) = A^VgW, EA{t) = X U) IJ$HRt = Xn Df i R t, 

3 = 1 

where 

^^(AW/A)^ (2.11) 
i=i 

is the mean of the mixture model of the durations of the different streams. 

Observe that we can write the cumulative inputs as linear drift plus compensated random Poisson fluc- 
tuation as follows: 

/OO pOO pOO (j) 
/ / rL t (s,u)N {ds,du,dr), (2.12) 
-oo JO JO 

/oo />oo />oo 
/ / rL t (s,u)N(ds,du,dr). (2.13) 
-oo JO JO 

After scaling time by T, we think of A$: = (A^(Tt),t > 0), j = 1, . . . ,p and A T : = (A(Tt),t > 0) for 
large T, as the cumulative inputs on large time scales. Thus, we consider a family of models indexed by the 
time scale parameter T and from now on we let the arrival intensities depend on T so that A^ := A^(T). 
If necessary, we let Xj(T) — >• oo as T — > oo (see (2.15)). Dependence of the arrival intensities on T means 

A, F D , b { £ and bjj as defined in (2.1)-(2.4) depend on T as well; however, notice that the tail indices of 
the distribution of the duration, namely a}y , remain independent of T. In practice, the fact that we focus 
on the stream at a particular time period, say [0,Tt], does not affect the tail index of the distribution of 
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the sessions duration, which is in accordance with our assumptions. For convenience, we often suppress the 
subscript T. 

Fix j, 1 < j < p and in the Tth model, let A^ B \t) be the centered and scaled cumulative input of the jth 
stream in [0,Tt], that is 

w ._ AU)(Tt)-\U)$v R Tt 

for a suitable aj(T) to be made precise below. Assuming lirriT->-oo X^TFfj (T) exists, the asymptotic 
behavior of A^J(t) as T — > oo, depends on whether the arrival rate is large, moderate, or small, relative to 
the tail of the duration. 

Theorem 2.1. (Mikosch et al, 2002; Kaj and Taqqu, 2008). 

For any 1 < j < p, consider the following three growth regimes of the arrival rate: 



lira X^TF^(T) = 



T-s-oo 



oo, fast-growth. 



Cj" , moderate- growth, (2-15) 
t 0, slow- growth, 

where Cj G (0, oo). (The form of the moderate- growth limit facilitates a simple expression of the corresponding 

limit process.) Assume that either F R G RV- aR ,a R > ctjy or E[(R^) 2 ] < oo. (If an < a^, the limit 
process is the same for all three growth regimes and the distinction among the growth regimes is irrelevant 
(Kaj and Taqqu, 2008, Theorem 4)-) 
(a) Under fast- growth, we distinguish two subcases: 
(i) IfE[(R^) 2 ]<oc, 

A(£(.)^E[(R^ 2 a^ u){1) B H u^ T^oo, 

where 

a«(T) = [X^T 3 F^(T)} 1/2 , 
<r? - 



b h u,W (a g)_i) (2 -aW)(3-ag)) : 
and B H (j) is a fractional Brownian motion with Hurst exponent 

= (3-a$)/2 G (1/2,1). 

(ii) If F R G RV- aR ,l < cq) <a R < 2, then 

A^(.)^Z a a )aR (.), T^oo, 

where 

a^{T)=Tb R {\TF { S\T)), 
r°° r°° r°° oo 

Z a/£, a Jt) = ] J o J o rL t(^ U ) N a ^\a R ( ds ^u), 



l/ag> 

- cos ^-g- ) 7, ^ (* D : / / L t (s,u)M aR>m (ds,du), 



na$\ 2r(2-a«) 



OO rOC 



and M aR ^ m (ds,du) is a a R -stable random measure with control measure 

m(ds,du) = ds ■ ap > u~ < - a n +1 ^du. 
Thus, the process Z q (t) is a R -stable and H^' -similar with 

HV> = (a R + 1 - a$)/a R G (l/a R , 1). 
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(b) Under moderate- growth 

fidi 



where 
and 



A^(.)^ Cj Y a ^(./ Cj ), T^oo, 
a (J)(T) = T, 



oo poo poo 



(c) Under slow- growth 



where 



Y a U)( t )= / / rL t {s,u)N a u) FR {ds,du,dr). 

' -oo JO JO 



4i ) (-)^^[(4 1) ) Q "] 1/ ^ , A Q(3) (-), T 



»«)(T) = b%\\^T), 



A (j) is an -stable Levy motion totally skewed to the right, which we can write as 



oo poo poo 



E[(R[ ') aD ] 1/a ° A Q ( 3) (t) = / / / url {0<s<t} N a U) tFR (ds,du,dr) 



oo JO JO 



ira 



oo poo 



l/ag> 

ag J (a$ - 1) j J-ooJo 



$\ 2T(2-ag>) 



-cos^- j - ) 



and M a u) m (ds,dr) is an a £ -stable random measure with control measure 

m[ds, dr) = dsFn(dr). 

Real network traffic consists of several distinct types and in this paper we are interested in the centered 
and scaled cumulative input of the superimposed streams in [0,T£], namely 

Acs[t) '= alT) ' ( 6) 

for a suitable a(T). In order to study the limit distribution of A cs (t) as T —> oo, let F, Ai, S be the subsets 
of indices of streams whose arrival intensities behave under the fast-, moderate-, and slow-growth regimes, 
respectively. 

Assuming that all indices belong to one of these three classes, consider the following scenarios. 

Scenario F: There is at least one stream whose arrival intensity satisfies fast-growth; i.e. F / 0. In 
this case, the aggregated stream's arrival intensity also satisfies fast-growth: 

XTF D (T) > \ {J) TF^ ] {T) -> oo, T ^ oo. (2.17) 

Scenario Ai: No stream's arrival intensity satisfies fast-growth, but at least one stream satisfies 
moderate-growth; i.e. F — and Ai ^ 0. Then, the aggregated stream's arrival intensity satis- 
fies moderate growth, since 

XTF D {T) -> c aD -\ T^oo, (2.18) 

where 



p 



OL D := f\ 0$ (2.19) 

3=1 

and 




(2.20) 
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Scenario S: All the stream's arrival intensities satisfy slow growth, that is S = {1, . . . ,p}. In this 
case, the aggregated stream's arrival intensity also satisfies slow-growth: 

\TF D {T) = \ (j) TF^(T) -> 0, T — > oo. (2.21) 

The different growth regimes in Theorem 2.1 are specified by the arrival intensity \W , and the distribution 
F^jp of the duration of the sessions of the jth stream. While A ^ = \V> (T) — >• oo as T — > oo, F^ does not 
vary with T . However, the growth regimes described in Scenarios T , M. and S are given in terms of the 
arrival rate A, and the distribution Fd of the duration of the sessions of the aggregated stream and here both 
A and Fd vary with T, as seen in (2.2). Therefore, we cannot directly apply Theorem 2.1 for the aggregated 
stream when 

A") /A = proportion of the sessions that belong to the jth stream, j = l,...,p, (2.22) 

are functions of T. Nevertheless, in the special case that these proportions are constant, Fd does not vary 
with T, and a direct application of Theorem 2.1 yields the following result. 

Corollary 2.2. Suppose that for all T (or at least for T large enough), the proportions X^'/X remain 
constant, j = 1, . . . ,p, so that 

where cud is given in (2.19). Let the Scenarios J- ', A4 and S take the place of the fast-, moderate- and 
slow-growth regimes. 

If F R G RV- aR ,aR > a D or E[{Ry) 2 ] < oo, then Theorem 2.1 holds for A cs (-), where d$ , F^y and cj 
are replaced by old, Fd and the constant c in (2.20), respectively. 

If Fr 6 RV- aR ,a.R < aD, the distinction among Scenarios T , M. and S is irrelevant, and limit results 
are discussed in Section J±. 

Implications of Corollary 2.2. This result provides a partial answer to the question of how much aggre- 
gation is required for traffic to be Gaussian at large time scales: Suppose that at least one traffic stream 
falls in the fast-growth regime, thus generating a cumulative input that can be approximated by fractional 
Brownian motion. When applicable, Corollary 2.2 implies that the superimposed traffic load can also be 
approximated by fractional Brownian motion. 

In the case that the traffic also contains streams that satisfy the slow- growth regime, Corollary 2.2 is 
somewhat counterintuitive due to the nature of the distribution tails of the two limit processes. Although 
these slow-growth streams produce cumulative inputs that are approximately stable Levy-motion when 
considered individually, with the inclusion of one single stream that behaves under the fast-growth regime, 
the cumulative aggregated input is approximately Gaussian. 

Moreover, a sufficient condition for the fast-growth regime of Scenario J- is that a single stream, say the 
jth one, satisfies fast-growth, even if all the other streams' arrival intensities do not follow a growth regime 
at all. In this sense, Scenario is a robust assumption. We will see that as long as one cvjj < an, the limit 
result of Corollary 2.2 is still valid. 

In real networks, there are arguably streams with large initiation rates. For instance, the arrival rates 
of http traffic must be large, since there are a large number of users constantly accessing websites and this 
translates into Scenario T . Furthermore, even though some studies report or assume session transmission 
rates have infinite variance, the assumption _E[(i?J 1 '') 2 ] < 00 may be justified by rate constraint mecha- 
nisms required for congestion control. Although assumptions always deserve rigorous scrutiny, Corollary 2.2 
provides a compelling explanation for the data example in Section 1. 

We now address more general assumptions which allow the conclusions of Corollary 2.2 to hold. While 
the assumption of constant proportions X^'/X may sometimes be reasonable, in general the proportions 
of sessions corresponding to the p independent streams are not constant over time. We may have that 
IhnX^/X exists or, more generally, that X^'/X £ (a, 6) C (0,1) varies with no limit whatsoever. Extending 
the conclusions of Corollary 2.2 to under weaker assumptions is the focus of the next section. 
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3. Behavior of cumulative load of aggregated streams 

We prove that the conclusion of Corollary 2.2 is still valid even when the proportion of the sessions 
corresponding to the p independent streams is not constant. Here is the result: 

Theorem 3.1. Assume that 

liminf \/ A (J VA>0. (3.1) 

3-a { D=ar> 

Then, the conclusions of Corollary 2.2 regarding the limit distribution of the cumulative input of the aggre- 
gated stream A cs {-) are still valid. 

Condition 3.1 implies that there exists d > such that for all T sufficiently large, there is at least one 
k = k{T) such that = arj and A( fc )/A > d. Roughly speaking, this means that the proportion of the 
traffic with the heaviest-tailed duration always remains greater than a positive quantity. 

All the limits in Theorem 3.1 follow from the convergence of the characteristic function of the finite- 
dimensional distributions (fidi chf) of the processes. Thus, let to > 1 represent the dimension, < ti, . . . , t m 
the times, and z\, . . . ,z m arbitrary real numbers; we need 

{m m 
*X^ ri *j( s ' w ) ( 1 - ZjrL tj (s,u), 
i=l J 3=1 

as defined in Proposition 5.6. 

From the second integral in (2.8), we can compute the partial derivative of L t (s,u) with respect to u, 
which yields 



g u (s-u,u,r) := ^-g\ = i I cxp I i y~] ZjrL tj (s - u, u) \ - 1 ) ^ z fc rl [0 . tfc] (s), (3.2) 

0U '<™> V j= i / fc =i 



where g u is the partial derivative of g(s,u,r) with respect to u. Moreover, putting together (5.4), (5.16), 
(5.18), the bounds in Lemmas 5.4 and 5.7, we get 



exp < i z 3 rL tj (s - u, u) > — 1 
i=i 



< 2^ \zj\Htj A«) ( r ( , < C < 1- (3.3) 



We will use three more relations in the proof of Theorem 3.1: For < r\ < 1, there exists a number 
T = T (77) > such that for T > T and 6 D (AT) > T , 

2»->{«-"V«"} > \^)I^ T \ ^ /T, (3.4) 

1 > ~ \F D (b D (\T)u)/F D (b D (\T)), u>T /b(XT), V ; 

< voT^-^u- 1 , u<T Q /T, (3.5) 

and 

¥^7^V7 < »nbn{Xrr°-^u-\ u < n/b D (\T). (3.6) 
F D (b D (\T)) 

We can readily derive (3.4) from Lemma 5.3. Both (3.5) and (3.6) follow from Markov's Inequality and, for 
example, Resnick (1987, Proposition 0.8) or Bingham ct al (1987, Proposition 1.3.6). 

Proof of Theorem 3.1. First, we will prove parts (b) and (c). For both parts, set 0<?7<a.D — 1 and 
< C < 1 such that 

, , . ia R , if F R G RV- aR , Ka R <2, 
°">+"< 1 + «U iiE [(R^n<oc. 
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Part (b). Under Scenario M, use a(T) = T and apply Proposition 5.6, yielding 

! m roc poo roo Fd(Tu) 

i'S~\z m A C8 (tj) > = jfj 9u(s - u,u,r)XTF D (T) - dsduF R (dr), 
~[ J-ooJo Jo D\ ) 

and if we can take the limit inside the integral as T — ¥ oo, performing afterwards an integration by parts in 
u gives 



oo roc poo 



g u (s - u,u,r)c aD ~ L u~ aD dsduF R (dr) 



-ooJO JO 

oo roo poo 



-oo JO JO 



g(s,u,r)EN°° (j) (ds,du,dr), (3.7) 



which is the log hdi chf of cY aD (-/c). Thus, it suffices to justify taking the limit inside the integral. 
First observe there exists a number To > such that for T > Tq, 

XTF D (T) < c 010 - 1 +n, (3.8) 

by the moderate-growth assumption. Together with (3.3)-(3.5) and a possibly larger To, the above implies 
that the integrand in the left side of (3.7) is bounded in {u > T /T} by 

m m 

B M , {>) (s,u,r) :=4(c Q -- 1 +7 ? ) U -«-( U -"V U ")EEl^l C l^l^ Au ) Crl+Cl [o^]( s )' 

j=l k=l 

and bounded in {u < Tq/T} by 

m m 

B M , (<) (s,u,r) :=2(c«- 1 +,)C- 1+ 'V I) ^^|z/|z t |^«-V^l [0 , ttl ( s )l (0 , 1) ( I1 ), 

3=1 k=l 

whenever T > Tq. Here we used the bound 

u c < (T /T) QD - 1+r 'u 1+c - QD -'', < u < Tq/T. (3.9) 
Therefore, (3.7) follows by the dominated convergence theorem, since for all T > Tq 



OO rOO fOC 



OO JO JO 



B M ,(»( s ,'u,r)dsduF R (dr) 

mm ( f 1 f°° 

< 4 ( cQD -i +?? ) £ ;[( i? (i))i+C]^^| Zj |C| Zfe | t J / u Z-<*n-n du + t i I u - aD+ ^du 

fe=l3=l Jl 

and 



/OO rOO rOO 
/ / B M ^(s,u,r)dsduF R (dr) 
-co JO JO 

mm „i 

<2( C «- 1 + )) )C- 1+ V^^|z J f|z t |t fc £p« 1) ) 1 «] / uC-^-'J 

3=1 k=l J ° 



du. 



which arc both finite by our choice of r\ and £. 

Part (c). Under Scenario <S, a(T) = bp (XT), so we use Lemma 5.1 and Mikosch et al (2002, Lemma 1) 
to get 

lim T/a(T) = oo. 

T-^oo 

Thus, it follows from the definition of L t (s,u) in (2.8) that 

lim L tT / a ( T )(sT/a(T) -u,u) = uh ot] (s). 



12 L. LOPEZ-OLIVEROS AND S.I. RESNICK 

Now, apply Proposition 5.6, perform the change of variables r i— y ra(T)/T, u M> uT/a(T), and use the 
scaling property in Lemma 5.4 to get 

fin 
3=1 



oo />oc poo 



oo J JO 



g u (s - ua(T)/T, ua(T)/T, rT / a{T))\a(T)F D (a(T)u)dsduF R (dr) 



/oo poo poo I I rn \ 

* I exp j z i rL t,Tla{T){sTla(T) - u,u) | - 1 j 



x ^ z k rl {0 j k) (s)\TF D (a(T)u)dsduF R (dr)i 

k=l 

and assuming we can take the limit inside the integral, the limit as T — > oo is 

/oo poo poo I [ jn \ m 

^ ^ i f exp r^z 3 url^ tj) (s) \ -1 j Y^Zkrl^t^u-^dsduFnidr), (3.10) 

which is the log fidi chf of E[(F£f > ) aD ] l l aD K aD (■). Therefore, we must justify passing the limit inside the 
integral. This is done as follows. 

First, by Lemma 5.2, there exists a number To > such that for T > To, 

XTF D (a(T)) < 2. (3.11) 

Hence, by taking a possibly larger To, (3.3), (3.4) and (3.6) imply that the integrand in the left side of 
(3.10) is bounded in {u > T /a(T)} by 

m m 

B Si(>) (s,u,r) :=8 M -« D ( U -"V U ")EEl z /NI^ Au)<r 1+ % itk] (s), 

j=l k=l 

and bounded in {u < To/a(T)} by 

m m 

B S , (<) ( W ) :=4Tr- 1+ VDEEl^l C l^l wC " aD "" ?4+Cl [o^]( s ) 1 (o-i)( u )' 

i=i fe=i 

whenever T > Tq and a(T) > To, using 

" C < (T /a(T)) aD ~ 1+r >u 1+i: ~ aD ~ v , 0<u< T /o(T). 

Therefore, (3.10) follows exactly as in part (b) from the dominated convergence theorem. 

Part (a). Under Scenario T and E[(n[ 1} ) 2 ] < oo, set a(T) = [\T 3 F D (T)]^ 2 . Use Proposition 5.6 and the 
change of variables r H> ra(T)/T, to write 

In E exp L^z m A cs (tj)\ = J ^ jf g u (s - u,u,rT/a(T))(a(T)/T) 2 p ?^ dsduF R (dr), (3.12) 

where 

(a(T)/T) 2 = \TF D (T) -> oo, T -> oo, 

by the fast-growth assumption. 
By (3.2), as T -> oo 



g u (s - u,u,rT/a(T))(a(T)/T) 2 = i\i ^ z j r ~J(fj L ti( s ~ u > u ) + (^)) ) E^^to.^W 



r 



HETEROGENEOUS TRAFFIC AT LARGE TIME SCALES 13 

Hence, assuming we can pass the limit inside the integral, we use Lemma 5.4 (hi) to write 



lim \nEexp{iy z m A cs (tj)} 

^— ^oo ^ — * 

3 = 1 

m m poc roo poo 

EE 2 ^ / / / r 2 L t] {s - u,u)l [Q , tk] (s)u~ aD dsduF R (dr) 

,_i J-oo JO JO 



3=1 k=l 



= -e[(r[ 1) ) 2 } 



i 



{a D - 1)(2 - a D )(3 - a D ) 

m J mm 



{ ill j m in 

EI>^- aD +E E *M$~ aD -^-h?- aD ) 
j=l k=l j=l fc=j+l J 

_. rn m 1 



j=i fe=i 

where the last line follows by rearranging of the terms in the sum, o^h(i) * s §i vcn m (5-10) and H = 
(3 — Q!£))/2. It remains to prove that we can take the limit inside the integral. 

Let < 7] < an — 1- We use (3.3)-(3.5) with £ = 1 and a possibly larger To, which imply that the 
integrand in Eq. 3.12 is bounded in {u > T /T} by 

m m 

j=i fc=i 

and bounded in {u < Tq/T} by 

^ >(<) := 2T q -- 1+ Vd E E |^*|« 1 - aD - , 'r 2 l [ o )tfc ](*)l ( o I i)(it), 

3=1 k=l 

whenever T > Tq. Here we used 

u < {T /T) aD ~ 1+r 'u 2 ^ aD - r i, < u < T JT. 
The result now follows by the dominated convergence theorem, since 

/oo />oo />oo 
/ / Bjr t ^(s,u,r)dsduF R (dr) 
-oo Jo Jo 

m m ( r 1 r°° l 

fCLEKi^VlEEl^l^ / u^^du + tjER^ u- a ^du\, 

j=l k=l lJ ° Jl J 

and 

/oo y* oo poo 
/ / Bjr {<) (s,u,r)dsduF R (dr) 
oo JO Jo 

m rn „\ 

< 2T a -- 1+ Vz) V V ^^[(i?^) 2 ] / u^-odu, 

3 = 1 k=l J ° 

and both bounds are finite by our choice of r\. 

Finally, still under Scenario J 7 , assume F R G i?V_ QB ,l < a R < 2. Set a(T) = Tb R (XTF D (T)). By 
Proposition 5.6, an integration by parts in u and the change of variables sh>s + u: 

»E Zm ^(^ 

3=1 
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= XTF D (T)F R (b R (XTF D (T))) / / g(s,u,r)ds - 

J-ooJo Jo ^b(r) F R {b R {XTF D (T))) 

= XTF D (T)F R {b R (XTF D (T))) {l { + /(r<e)}, (3.13) 

where e > and we split the integral into three parts according to the domains of integration {u > e, r > e}, 
{u < e,r > e} and {r < e}, respectively. To establish the limit result, we will take lim e _>o limx^oo on both 
sides of (3.13). 

Fix e > and start with the first integral. Let 

, , , s ( F D (Tdu)\ ( F R (b R (XTF D (T))drY 
Vxiflu, ar) — 

v(du, dr) 

G(u,r) 



F D (T) ) V F R {b R (XTF D (T))) 
anu~ aD ' du a R r~ aR dr, 



1 f°° 

— / g(s,u,r)ds, 
ur J_ 00 



which allows writing 



oo poo 



I(u>e,r>e) = / G{u, r)v T {du, dr) . 



The fast-growth regime, regular variation of F R , (5.4) and Billingslcy (1999, Theorem 2.8) imply vt — > v 
as T —> oo. Moreover, G(u,r) is jointly continuous and it follows from Lemmas 5.4 and 5.7 that \G(u, r)\ < 
dp Y]j— 1 \zj\tj < oo, where do is a positive constant. Therefore: 



T-s-oo 



OO r OO 



Jim J (u>£ , r>e ) = / / G{u 1 r)v{du,dr) 



OO />OQ z 1 OO 



g(s, u, r)ds ■ a D u- {aD+l ^du ■ a R r- {aR+l ^dr (3.14) 



Now let < C, ?y < 1 such that old + f] < 1 + C < By Lemmas 5.4 and 5.7, there exists > such 
that 



\I(u<e,r>e) \ — 



J 

Furthermore, by fast-growth and regular variation of F R 



J ST\ \i r i +C FR(P R (XTF D (T))dr) f f 1+c F D (Tdu 
'ff "./ r F R (b R (XTF D (T))) L U F D {T) 



Je F R {b R {XTF D {T))) J t chr-I-C 

Similarly, integration by parts and (3.5) with T >Tq such that e < Tq/T yields 

F D (T) 1 S Vo ^d(T) 



< (1 + C)^T " D - 1+ " / u^^du 

^ (i + 0^r;°- 1+ " 

~ 1 + £ - a/? - rj 
where we used the bound (3.9). Thus 

limsup | J(„ <e , r > e ) | < constant ■ e 2{1+ °- aR - aD ~^ , (3.15) 

T->oo 

where the exponent of e is positive if we additionally let (a R + ajj + n)/2 < 1 + C- 
For h r <e)i use again Lemmas 5.4 and 5.7 to get a d\ > such that 

\T \*A V^l |4 f 2 F R {b R {XTF D {T))dr) [°° Fp (Tdu) 

l (r<£) ' SS |ZjZfc ' ^o r F R ( bR (XTF D (T))) I " ~ Fb^T) - ' 
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Analogously to the bound for I( u <e,r>e)i it can be readily shown that there exists Tq > such that for T > Tq 

r2 Fn(b R (XTFp(T))dr) < 3 m T«« ^ 2 _ aR 



F R (b R (XTF D (T))) ~ 2-a R 



and 



whence 



^(Trfu) f 1 2 F D (Tdu) [°° F D (Tdu) 3» D T£° a D 
u(t k Au)—=—— < u - + t k / u —=—— < — \-t k - 







lim sup I I(r<e) I 5= constant ■ e 2 Qi? . (3.16) 

T^oo 

Also, by fast growth and the regular variation of F R 

\TF D (T)F R (b R (\TF D {T))) -+ 1, T oo. (3.17) 
Finally, we can put together (3.13)-(3.17) to write: 



lim lnSexp < ZjA cs (tj) > = lim lim ln_Bexp < i ZjA cs (tj 

{ 3=1 

) 

flr(s, w, r)ds • a D u- (aD+1) du ■ a R r- {aR+l) dr. 



e->OT^oo 

3=1 3 = 1 



oo /* oo /-oo 



oo JO JO 



□ 



4. Remaining choices of old and a R 

We first study what happens if a R < Od, namely, a R < oqj for j = 1, . . . ,p. Consider the centered and 
scaled cumulative input of any of the streams, say the first one, for simplicity. Set the normalizing term to 
q(i)(T) = b R (XT) (and the same for all other streams). Write 

Y POO />00 />00 f]\ 

A -(*) = t-7T7^ / / / rL Tt (s,u)N (ds,du,dr) 



OO /» OO /"OO 



o (1) 

rL>Tt(Ts,u)N (Tds,du,b R {\T)dr). 



'-oo JO JO 

First, note from the definition of Lt(s,u) in (2.8) that 

lim L T t(Ts,u) = ulr tl (s). 

T^oo 

Thus, analogous to the proof of Proposition 5.6, it can be shown that the log fidi chf of A^J is 

A« - 
= —j—\TF R (b R (\T)) 

x Z^/ 00 / 00 * ^cxpjzrfjz.L^^,^! -lj p^ k L Tt] (T S ,u)d S F { D 1} (du)^M^dr. (4.1) 
Now observe that, provided we can take the limit inside the integral 
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/oc poo poo I m \ m 

J J i I exp | ir^Zj-wlfo^s) \ 1 I ^z k ul [0:t]] (s)dsF£ } (d 



u)r- aR dr. (4.2) 



Since | A (1 > /A| < 1, then 



AW 
A 



/OO />oo /»oo / m \ m 

y J i I exp j ir^z J ul [0! i 3 ](s) > 1 j ^ z k ul [0 ^ t:j] (s)dsF^\du)r- 



l dr 



->-0, T^oo, (4.3) 
which yields the following result. 
Theorem 4.1. Lei * := *t &e ifte /id? eft/ of 
p \U) 

52—E[(DP)''*] 1 ' a *A aR (t) 

= ((— ™) ^^y^cr^-^-^ {iA> 

where for each T, A QR (-) is an a R — stable Levy motion totally skewed to the right with index a R and 
M aR (ds,du) is a R — stable with control measure m(ds,du) = dsFjj(du). Then, 



lim Iln^exp ji^Zj-AWfo) j -ln* tlj .. 



..., tm (z u ...,z m ) } =0. (4.5) 



In addition, if for j = 1, . . . ,p, the limits 



w {3) := lim (4.6) 

T— s-oo 

exist, then the fidi chf of A cs (-) converges to the fidi chf of the process defined by (4.4), with w^i and 
T,%i wU)f d ] (-) replacing \W/X and F D (-). 

Proof. By the independence of N^\j = 1, . . . ,p, 



In E exp < i ZjA cs (tj 



i=i 



/OO pOO pOO I m \ m 

/ / i exp I iry2zjul[ , tj ](s) > 1 V] z fe Ml [0 . tj . ] (s)dsi^(cZ'u)r~ QR c/> 
00 ^ Jo \ [ j=i J / fc=i 

p 

= ^ e C*)(T)-+0 T^oo. (4.7) 
i=i 

Analogously to the proof of (5.13), the second integral in (4.7) is equal to hx^fti,...,t m (^ij ■ ■ ■ , z m)- Thus, it 
only remains to justify taking the limit (4.3). 

Let < C < C < 1 an d < 77 < 1 such that 1 + ( < a R ■- n < a R + 77 < 1 + (' < a D . Similarly to (3.4) 
and (3.6), there exists T := T (r]) > such that for T > T and b R (XT) > T , 

F R (b R {XT)r) < f 2r-°« {r~* V r"} , r > T /b R (XT), 
F R (b R (XT)) -[uRbHiXT) *- 1 ^- 1 , reR. 
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Together with Lemmas 5.4 and 5.7, this implies that the integrand in the left side of (4.2) is bounded in 
{r>l} 

m m 

B(» :=2 1 -V + V- a « + ^^|z J f|^|l [0jtfe] ( S )l [1 ,oc)M, 

j=i fc=i 

and bounded in {r < 1} by 

m m 

S (<) :=2 1 ^V^- 1+ V^'/'-^-^^|z J f'|z fc |l [0 , tfc] ( S )l ( o,i)M, 

whenever bn(XT) > T . Here we used 

r c ' < {T /{b R (XT)) aR - 1+,1 r 1+<; '-"D-fl . 

By our choice of £, £' and 77, both bounds are integrable and we can use dominated convergence to prove 
the result. □ 

In principle, it also is possible to have old = olr. However, we cannot say much except in the special 
case old = ctR = 2, in which case the limit process is a Brownian motion provided (4.6) holds. We refer the 
reader to Kaj and Taqqu (2008, Theorem 4) for the formal statement of this case. 

5. Technical proofs 

This section contains a collection of technical results needed for our proofs. The first lemma establishes 
bounds for 6d(-) = (1/Fd)*~(-) which yield bo(XT) — > 00. This is not immediate since the function bo 
depends on T. 

Lemma 5.1. The quantile functions given (2.3) and (2.4) satisfy the following inequality. 

(/ b%\ P X^T) > bp (XT) > \) b$(\®T), T > 0. (5.1) 
3=1 i=i 

Hence 

bjj(XT) — > 00, T — > 00. 
Proof. Since Ffj is decreasing for all j, then 

Fd ( V&gV«T)| <^(A«/A)4 j) (^(pA«T)) 
\i=l / 3=1 

< f^(X^ /X)(pX^T)~ l 

3=1 

= (XT)-\ 

Thus, the left side of (5.1) follows. 

On the other hand, since Fd is right continuous, we have for each j = 1, . . . ,p: 

(X^/X)F^(b D (XT)) < ^(AW/A)4 fe) (6 D (AT)) < (AT) -1 , 

k=l 

whence 

F%\b D (XT)) < (X^T)- 1 . 

Therefore, the right side of (5.1) follows. □ 

The distribution F D = ~£ p j=1 (A^ /X)F^ of session durations of superimposed streams is a function of T 
since A^ and A depend on T. Nevertheless, Fd behaves as a regularly varying function. 
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Lemma 5.2. Under the assumption (3.1) 

lim lim \TF D (b D (\T)u)=u- a °, u>0, (5.2) 

and therefore, in Af + (0, oo], 

£|™ A c^-^W, (5 . 3) 

Proof. Note that Fd G RV- an for each fixed T. However, because Fd varies with T. the limit is not 
straightforward. 

Fix an arbitrary u > 0. We start by writing 



F D (Tu) _ Ej^W^'AlF^Tu) 



=:5 + ^ (A«/A)Q, 



and additionally write 



E, :a g) =QD (A0')/A)F«(T u ) E fa ^(A*/A)F«(Tt 
Bi+ E ( AW /A)i?2,. 



Thus, the first limit in (5.2) will follow by proving B\ — > u a ° , i?2j — > and Cj — > as T — > oo, for all j 

(i) 

such that > an. 

First, by Potter's bounds applied to the regular variation of each Ffj , we have as T — > oo, 

£ rQ o>- QD (A W /A)^(T) 

E Jiag)=ai ,(AW/A)^(r)«-« 
( \ (j) (j) 

Now, consider B2J for a}y > ap. Choose an arbitrarily large z > mm.{u~ a D ,u a o }. By regular 
variation, for T sufficiently large: 

Fg\Tu) > _ F D k \T) ^ 
Fj»(Tu) ' Fj?(T) ' 



(k) 

for all k such that a D — ao- In addition 

^(Tu) -a«> -1 P d( T ) a«> -1 

-^7-7 ->u -z , > u d —z. 

F%\T) F<£\Tu) 

Furthermore, the assumption (3.1) means there exists d > such that for all T sufficiently large, there is 
some k' := k'(T) such that ' = a D and A^ fe '/A > d. Hence for T sufficiently large: 

This shows that B^j can be made arbitrarily large for T sufficiently large, whence B2J — > as T — y 00. 
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Similarly, consider Cj, and 



>dz(u aD -z- L ). 



This shows that Cj can be made arbitrarily large for T sufficiently large, which completes the first part of 
the Lemma. 

For the second limit in (5.2), recall that z < bo{XT) iff 1/Fd(z) < XT for each T. For e > 0, setting 
z = b D (XT)(l - e) and z = b D (XT)(l + e) yields 

F D (b D (XT)(l + e)) 1 F D {b D {XT){\ - e)) 

F D (b D (XT)) ~ XTF D (b D (XT)) ~ F D (b D {XT)) 

Letting T — > oo and using Lemma 5.1 and the first limit gives 

(1 + e)~ aD < = — < (1 - e)~ aD . 

[ ' ~ XTF D (b D (XT)) ~ [ 1 



Because e is arbitrary, then 



lim XTF D {b D (XT)) = 1. 

T— ¥00 



Therefore 



lim XTF D {b D {XT)u) = lim XTF D (b D (XT)) lim F bM^M = u -«n , 



The final statement about vague convergence follows the proof of Resnick (2007, Theorem 3.6). □ 
Even though Fd depends on T, a version of Potter's bounds holds. 

Lemma 5.3. Let S > 0. Under the assumption (3.1), there exists Tq = Tq(S) > such that for all T > Tq, 
Tu > T Q : 

§^<(l + ^-^max{ U -V}. 

Proof. Observe 



V 



Fd{T) E^i(A^')/A)4 j) (T) 3 l, F«(T) ' 

By Potter bounds (See e.g. Bingham et al, 1987, Theorem 1.5.6), for all j = 1, . . . , p there exists Tj = Tj(S) 
such that 

Fd , < (1 + *)u -a ° maxlu- 5 , u s } < (1 + 5)u- aD max^" 5 , u 5 }, 

F«(T) " 

for T>T j7 Tu > Tj. Therefore, the result holds for T = \J P ]=1 Tj. □ 

We now study L t (s,u), as defined in (2.8). 

Lemma 5.4. The length of the subinterval of [0,t] during which the session (s,u, r) transmits data, namely 
Lt{s,u) in (2.8), satisfies the following properties: 

(i) Scaling property: For C > 0, 

CL t (s,u) = L C t(Cs,Cu). 

(ii) Bounds: 

L t (s,u) <tAu. (5.4) 
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(iii) Integrals: For 1 < 7 < 2 and nonnegative t±,t2, 

/oo 
L tl (s, u)ds = uti, 
-00 

and 

/>O0 />oo 

L tl (s — u,u)l[ 0tt2 ](s)u~ J duds 



00 Jo 



1 



( 7 -l)(2- 7 )(3-7) 
Proof. The scaling property and the bounds follow directly from (2.8). 

Now, the first part of Property (iii) is readily checked by using the first integral in (2.8) after reversing 
the order of integration. Finally, the second part of Property (iii) can be derived by writing 

0, s < or s > u + tx, 

s, < s < u Aii, 

tl, t\ < s < u 1 

u, u < s < t\, 

ti — s + u, tiV ti < s < u + ti, 



L tl (s -u,u) = < 



and integrating accordingly. Observe that the four regions in which L tl (s — u, u) is nonzero correspond to 
those of the basic decomposition in Mikosch et al (2002, Eq. 4.1). 

□ 

The next lemma helps obtain approximations to the cumulative input of the aggregated streams. 
Lemma 5.5. For any a,T > 0, we have 

III o 

rLxt(Ts, u)N(Tds, du, dr) 



-{A{Tt)-X i i Dl i R Tt) = - 
a a 



-00 jo Jo 

00 />oo />oo 



-00 JO JO 
00 />oo poo 



00 Jo Jo 



rL t (s, u)N(Tds, Tdu, dr) 
rL t (s, u)N(Tds, Tdu, (a/T)dr). 



Proof. All the relations here follow from several ways to change variables in (2.13) and using the scaling 
property of L t (s, u). See Lemma 5.4. □ 

Our limit theorems are proved by verifying convergence of finite dimensional distributions for various 
processes. The following is required. 

Proposition 5.6. For arbitrary m > 1, < ti, . . . , t m , and real z\, . . . , z m , define 



fJ 



(s,u,r) = exp< i^z.jrL t] (s,u) > - 1 - ZjrL tj (s, u), 



(5.5) 



and 



h(s,u,r) = i l^cxp z i url (o,ti)( s ) j l j Zkl lo,t k ] (s)ru a ° . (5.6) 

(a) For any a,T > ; the characteristic function of the finite- dimensional distributions (fidi chf) of the 
process {(l/a)(A(Tt) — Xfir)fJ.RTt);t > 0} is given by 



In E exp < i Zj 



~(A(Ttj) - Xfi D fJ>R.Ttj] 
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g(s, u, r)EN{Tds, Tdu, {a/T)dr) (5.7) 
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(•OO />OG pOO 

g u (s - u, u, r)\TF D (Tu)dsduF R {{a/T)dr). (5.8) 

-ooJn Jo 

where g u is the partial derivative of g with respect to u. 
(b) The fidi chf of the limit processes in Corollary 2.2 are given as follows. 

(i) The fidi chf of the limit process under Scenario J- and E[(T6p) 2 ] < oo is given by 

!m 

.. m m 

= ^[(^lEE (*r + *f - \U -*i| 2 ) , (5.9) 

j=l fe=l 

where Bh is fractional Brownian motion with 

2 _ 2 

~ (a D -l)(2-a D )(3-a D y (5 - 10) 

and H = (3 — ar>)/2. 

(ii) The fidi chf of the limit process under Scenario T and Fr £ RVL aR , 1 < < 2, is given by 

oo poo poo 



g(s,u,r)EN™, aR (ds,du,dr). ( 5 - U ) 
(Hi) The fidi chf of the limit process under Scenario M. is given by 

! m pOO pOO pOO 

iJ^ Zj cY aD (tj/c) \ ^c^" 1 J J J g(s,u,r)EN™ tFR (ds,du,dr). (5.12) 

(iv) Finally, the fidi chf of the limit process under Scenario S is given by 



{ rn pOO pOO poo 

i^E^r^K^) } = J oo j Q I 



h(s,u,r)dsduFfj(dr). (5.13) 



Proof. Given (5.7), (5.8) is readily derived using integration by parts and the change of variables sh-s + u. 
Moreover, (5.9) follows from the fact that Bh is fractional Brownian motion. The remaining parts are a 
consequence of the following property of Poisson random measures (See e.g. Rosihski and Rajput, 1989): 

InSexp ji J f(x)£{dx)\ = J (e if & - 1 - if{xjj E£,{dx), 

if 

r (f 2 (x)A\f(x)\)E£(dx) <oo. 



For now, let us focus on (5.7), (5.11) and (5.12). By Lemma 5.5, the exponent in the left side of (5.7) 
and (5.11) is of the form 



/oo />oo />oo nl 
/ / ^2zjrL tj (s,u)l(ds,du,dr), (5.14) 
-oo JO JO - =1 

for a PRM £, while the exponent in the left side of (5.12) is c QD_1 times (5.14), using Lemma 5.4 and the 
change of variables s >->• s/c,u i-> u/c. Thus, it suffices to check that 

2 



/oo />oo />oo / m \ 



E Z k rL t k 0, i 



fc=l 



E£(ds, du, dr) < oo. (5.15) 



Bounds and integral results for Lt(s,u) in Lemma 5.4 (ii) and (hi) needed. 
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First observe that 



oo poo poo 



oo J0 J0 
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EN(Tds,Tdu, (a/T)dr) < — \zj\X/iD^Rtj, 



which proves (5.7). 

In order to prove (5.11), split the corresponding integral (5.15) into two parts /(<) and -?(>), according to 
the two domains of integration £>(<) = {ur < 1} and = {ur > 1}. This yields 

oo poo p^-/u 
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2 — a.R V aft - a D 1 - a R + a D J ' 



and 
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r 2 L tj . (s, ujLt, (s, u)EN^ aR (ds, du, dr) 
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whence (5.11) holds. 

Similarily, split the integral (5.15) corresponding to the process (5.12) into two parts J(<\ and J(>), 
according to the two domains of integration -D(<) and £*(>), which yields 

jn m „oo poo pr^ 1 

J«) < EE l z J' Zfc l / / / r2L t 3 ( s ^ u ) L t k ( s ^ u ) EN a D ,F H ( d s,du,dr) 
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and 
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J = 1 J-ooJO Jr- 1 
m | | 

<rJw£B[ (fl (i))« D ]. 

r— f old — 1 



This proves (5.12). 

Finally, the exponent in the left side of (5.13) is 



,oo ,oo .oo m o ^ 

* / / / E^' ur V%]( s )^ctt,,F H (dMMr)- 
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Analogous to the proof of (5.12), it is readily shown that 
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EN™ D FR (ds,du,dr) < oo, 
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whence the left side of (5.13) is equal to 



/oo />oo poo I f 171 ^ \ 

Jo Jo l exp {*^ ZjUrl [°^]^' ) J 1 " IZ z J Mrl [o.t J ]( s ) j EN a D ,F R (ds,du,dr). 

The result now follows after an integration by parts in the variable u. □ 

Finally, the following result is used to get upper bounds for some integrands throughout the proof of 
Theorem 3.1. 



Lemma 5.7. For < C < 1 cmd x <E 



\e lx - 1| < 2 1 - < |x| c , (5.16) 
\e ix - l-ix\ < d ( \x\ c+ \ (5.17) 



where d^ > 0, and for real numbers x±, . . . , x m : 

n \ £ m 

£n <En c - ( 5 - 18 ) 

•3=1 ' j=l 

Proof. Without loss of generality, fix x ^ 0. Define / : [0, 1] -)■ K, /(C) = (1 - C) In 2 + ( In \x\ . We can 
readily check that 

ln|e M "l|</(C), C - 0,1, 
by taking logarithms in both sides of \e tx — 1| < 2 A |x|. Since /(£) is linear in £, /(C) is either nondecreasing 
or nonincreasing on [0, 1]. Hence, (5.16) holds. Using a similar strategy, we can prove (5.17). 

For (5.18), assume without loss of generality that < \xi\ < ja^l, thus < |a;i/x2| < 1. By Bernoulli's 
inequality (see e.g. Mitrinovic and Vasic, 1970, p. 36): 

(1 + \xi/x 2 \) C < 1 + C\xi/x 2 \ < 1 + |xi/x 2 | c . 

Multiplying both sides by \x 2 \^ proves (5.18) for m = 2 and the proof for general m follows by induction. 

□ 
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